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Abstract. A model with Holstein-like electron-phonon coupling is studied in the 
limit of adiabatic phonons. The phonon distribution is anharmonic with two degenerate 
maxima. This model can be related to fcrmions in a correlated binary alloy and 
describes microscopic phase separation. We discuss the weak and strong electron- 
phonon coupling limit and present a qualitative phase diagram. In terms of the 
phononic displacements it consists of a homogeneous, an alternating, and a disordered 
phase. There is a first order phase transition between the homogeneous and the 
alternating phase, and second order phase transition between the alternating and the 
disordered phase. The opening of a gap inside the disordered phase is treated by a 
dynamical mean-field theory. 

PACS numbers: 71.30.+h, 71.38.-k 
1. Introduction 

Adiabatic (or static) phonons have been discussed in the literature by many authors 
[1112101 HI- They can be considered as a first step towards a full treatment of the electron- 
phonon interaction in a many-body system. The latter is plagued by a number of 
difficulties like the unrestricted number of phonons at each lattice site. These problems 
are partially avoided by the adiabatic limit, where the phonons are described by classical 
degrees of freedom. 

In general, the electron-phonon interaction can be successfully treated in a self- 
consistent perturbation theory, as long as the corresponding coupling is small. This 
approach is known as the Migdal-Eliashberg theory E] • At strong coupling, however, 
this theory experiences problems known as the breakdown of the Migdal-Eliashberg 
theory. In the case of the Holstein model it was argued by Benedetti and Zeyher that 
this breakdown at large electron-phonon coupling is associated with the appearence of 
an effective double- well potential for the atomic positions [7J. It is conceivable that 
the appearence of a degeneracy in the atomic positions is a source of problems in a 
perturbative approach like that of Migdal and Eliashberg. 

The aim of this work is to discuss the physics of anharmonic adiabatic phonons 
with two degenerate maxima of the distribution at each lattice site in terms of a simple 
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model. We ignore the electronic spin which is of no interest in our study and concentrate 
only on the scattering of spinless fermions. Our discussion includes the derivation of 
an effective Ising-spin model, where the spins represent the degenerate maxima of the 
distribution, and the application of the dynamical mean-field theory, which corresponds 
with the infinite-dimensional limit of the model. Our anharmonic model can also be 
considered as the strong-coupling regime of the adiabatic Holstein model. 

The paper is organized as follows: The model is discussed in section 2 and an 
effective distribution of the adiabatic phonons is derived in section 3. In the limit 
of discrete phonon degrees of freedom (correlated binary alloy) this is related to a 
distribution of Ising spins (section 3.1). The latter is studied in the weak- (section 4.1) 
and strong-coupling (section 4.2) limits. In section 5 we apply the dynamical mean- 
field approach to our model. Finally, the results of the two approaches are discussed in 
section 6. 



2. The model for adiabatic phonons 

The Holstein-like model with spinless fermions coupled to adiabatic anharmonic phonons 
is defined by the Hamiltonian iff and the thermal phonon distribution Pq at inverse 
temperature (3: 

H{ = C t C r' ~ H (/■* ~ 9 X r)4 C r > ^o(^r) OC exp 

<r,r') r 

cj (c r ) are the fermionic creation (annihilation) operators and the real variable x r 
represents the phonon degrees of freedom. The latter describes the displacements of 
the atoms at sites r, assuming that it is thermally distributed according to P . The 
fermions feel the displaced atom as a one-body potential —gx r . For harmonic phonons 
we have a Gaussian distribution 

Ph(x r ) oc exp (-/fa;!) > ( 2 ) 

instead of Pq, with x r = as the position with maximal weight. The anharmonic 
distribution P , on the other hand, has two degenerate positions with maximal weight, 
namely x r = ±1. The coupling to the fermions, however, can break this degeneracy, 
leading to homogeneous and inhomogeneous equilibrium distributions of the atomic 
positions on the lattice. 

Although the model defined in equation looks like a model for Anderson 
localization this is not a correct interpretation. The reason is that a model for Anderson 
localization requires the averaging of the Green's function with respect to Pq. This is 
not the case for the adiabatic phonons: A grand-canonical ensemble of spinless fermions, 
coupled to adiabatic phonons, is defined by the partition function 

Z = Tr e~ pH = f Tr f e~^ Hi H P (x r )dx r (3) 



(i) 
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which represents an annealed average of the fermionic system. The fermionic Green's 
function then reads in the temperature formalism 



G 



r,t;r',0 



-Tr 



e-^ H c r e- tH cl 



1 

Z 



Tr f 



e -09-t)H f£ . e -tH f J 



C r'\ n^W^r 



(4) 



The anharmonicity is related to previous studies where an effective anharmonic 
(double- well) potential for the electrons was found jjj- Here the main idea is that the 
tunneling of the fermions is much faster than the motion of the atoms between the two 
maxima of the distribution. Therefore, the dynamics of the phonons is negligible and 
adiabatic phonons can serve as phononic degrees of freedom. 



3. Effective Phonon Distribution 

The trace with respect to the fermions can be performed in equations (J3J) and since 
the fermions do not interact directly in the Hamiltonian iff jSj. It gives the determinant 
of a N x N matrix, averaged with respect to the distribution P , such that the partition 
function reads jH] 

"det [l + e ^ +i - 9x) ] l[P (x r )dx r . 

r 

N is the number of lattice sites and t is the hopping matrix. The integrand of Z is a 
positive expression and can be considered as an effective distribution of the adiabatic 
phonons that replaces P : 

1 



P({ Xr }) 



-det 



1 + e 



f3(fi+t-gx) 



(5) 



3.1. Ising-Spin Representation 

For the evaluation of the effective phonon distribution given in equation (jSJ) we consider 
the limiting case U ^> 1. Then the distribution Pq selects 



S r = ±1, 



where S r is formally an Ising spin which describes discrete displacements of the atoms. 
This case can be considered as a correlated binary alloy (CBA), where the correlations 
are mediated by the fermions through the determinant in equation @. The partititon 
function is now a sum with respect to Ising spins: 

Z = Yl z (i S r}) with Z{{S r }) = det [l + e ^ +t - gS) ' . 

{S r =±l} 

Then the fermionic Green's function is a resolvent, averaged with respect to the effective 
phonon distribution 0: 

(iuj + H + i- gS) \ 

/ Ising 
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where (...)i S mg is the CBA distribution 

P({S r }) 



det 


1 _|_ e 0(n+i-gS) 




£{s r } d et 


1 _|_ e /3(M+t- 


gsj 



(6) 



It should be noticed that the CBA distribution is not Z<i invariant, i.e. not invariant 
under a global flip of the spins or the coupling constant g, in contrast to the distribution 
Po- 

4. Approximations of the CBA Distribution 

The effective phonon distribution given in equation (JHJ) can describe different types 
of order of the displacements of the atoms, depending on the inverse temperature f3. 
At high temperatures we expect a paramagnetic (disordered) distribution, and at low 
temperatures some kind of order. Some insight can be obtained by evaluating the 
distribution in the asymptotic regimes of weak as well as strong coupling. More details 
of the calculations are given in Appendix A. For further simplification of the results we 
restrict our interest to the low-temperature regime, i.e. to (3 ~ oo. 

4-1. Weak- Coupling Limit 

The weak-coupling limit is obtained from an expansion for small values of g. According 
to Appendix A the regime describes an uncorrelated binary alloy with the distribution 



p w ({Sr}) = n 



-/3ghS r 



with h 



Adv. 

©0" + ^))7^ (7) 



where e(k) is the dispersion of the tunneling term t. Thus the fermions create a 
homogeneous magnetic field gh for the Ising spins with < h < 1. This favors the 
atomic position x r = —1: 

(S) = -tanh (flgh). 

For a low fermion density (i.e. for // + e(k) < 0) there is no magnetic field (h = 0). 
Then the Ising spins are paramagnetic (disordered without correlations). 

4-2. Strong- Coupling Limit 

The strong-coupling regime enables us to apply an expansion in terms of the fermion 
hopping. As shown in Appendix A we can write for the effective phonon distribution 

exp{^ i [E 1 (S r + S r/ )+E 2 S r S r/ }} 



P s oc n [l + e^S r ) 



with coefficients 



*Pv-Pg 



^Pfi+Pg 



'I + e Pv-p g y (i _|_ e Pv+p g y 



'P 2 



^Pv-Pg 



^P^+Pg 



+ 



1 _|_ e Pv-Pg)2 (i _|_ e Pn+Pg)2 



and 



e w smh{j3g) 



Agl + e 2 ^ + 2efr* cosh(j3g) j ' 
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This gives in the low-temperature regime {(3 ~ oo) 



E x ~ 0, £ 2 



for -5 < n < g 
otherwise 



and for the distribution density 

r i 

exp { - /3[£(ji + g)S r /2 + P £ S r S rl /Ag] } 



P, ~ < 



(r,r') 



exp ( - /3<7 £ 5r 



for n < —g 
for -g < /!< g 

ioi g < ji . 



(8) 



In terms of the Ising spins there is a paramagnetic phase for /x < — g. For — g < /j, < g the 
competition of the magnetic field term and the antiferromagnetic spin-spin interaction 
leads to a first order phase transition between ferro- and antiferromagnetic states, at 
least for low temperatures. A simple mean-field approximation reveals that the Ising 
groundstate is an antiferromagnet for /i < dP/g — g and a ferromagnet for // > dP/g — g. 
Finally, for /i > g there is always a ferromagnetic state with (S) < 0. These three phases 
are shown in the low-temperature \x — g phase diagram in figure H 



5. The DMFT equations 

The spectral density of the itinerant fermions can be evaluated conveniently by the 
Dynamical Mean-Field Theory (DMFT), originally developed for the Hubbard model 
[31 EIJ \H\ ■ The main idea of the DMFT is the infinite-dimensional limit of the system. 
This limit is characterized by a simplified scattering process in comparison with the 
situation in a finite-dimensional system: the scattering of the same pair of fermions is 
very unlikely due to the large phase space. In the case of our adiabatic phonon system 
it is very unlikely that a fermion scatters with the same phonon more than once. This 
enables us to consider only the effective scattering on a single site. A consequence is 
that the effective correlation between the phonons, discussed in the weak- and strong- 
coupling regime, is excluded in the DMFT [T5j . 
Formally, the Green's function reads 

G.= f : * (9) 

with the Matsubara frequency u n and the density of states p of free fermions on the 
lattice. Moreover, the Green's function can also be written as j3] 

G » = r=i + + r-.'rr • ( 10 > 

Equation (fTUj) yields the self-consistent equation of the self energy S n : 

1 1 



^ = ± - V^ 2 + V + 4^(1 - 2w )G'\ (11) 
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where the coefficient u> depends on the self energy and G n as 




From these three equations the self energy S n and the Green's function G n can be 
determined numerically [3]. 

For a special form of the density of states p the integral in equation (jUJ) can 
be evaluated explicitly. This simplifies the numerical effort for solving the equations 
substantially. An example is the case of a semi-circular density of states p(e) = 
y/4:t 2 - e/2-nt 2 which gives fTJ] 



with ( n = iuj n + p — E„. wo is determined numerically using equations (fTT| and ((T2j). 
Substituting E n = iu n + fi — G^ 1 — PG n in equation (fTU|) yields a cubic equation for the 
Green's function |14j : 



with i] = iu n + p. Analytic continuation {iu n — > a;-|-iO + , G n — > G(u>)) and the condition 
Im G(u) < leads to the fermionic density of states D{uj) = — Im G(u))/n. Our results 
are shown in figure 121 

6. Discussion of the Results 

The phase diagram of our model is complex and has been studied by two different 
mean-field approaches, a classical one for the Ising-spin representation (i.e. the CBA) 
and the DMFT. It consists of para-, ferro- and antiferromagnetic phases in terms of 
the Ising spins (cf. equation ©)• The latter correspond to displacement configurations 
of the adiabatic phonons: disordered phonons (or phonon liquid), homogeneous and 
alternating (charge-ordered) phases. There are phase transitions between these phases. 
According to our mean-field calculations, a first-order transition appears between the 
metallic homogeneous and the insulating alternating phase. The transition between the 
disordered and insulating alternating phase is second order. The disordered phase is 
most likely also insulating due to Anderson localization of the fermions. The weak- 
coupling behavior of equation indicates a disordered phase for \i < —2di and a 
homogeneous phase for \i > —2di. These results are summarized in the qualitative 
phase diagram of figure HJ 

Our DMFT results do not agree very well with the picture obtained from the CBA. 
This might be a consequence of the fact that the DMFT is a local approach, at least 
in the version discussed in this paper |13j . There is a critical coupling g c that indicates 
the opening of a gap for g > g c . In our model g c changes with \i and T (cf. figure EJ- 
Large positive values for /i imply g c rs 0.5 for hopping rate t = 0.5. Decreasing /z results 
in a reduced value of g c . The lower the temperature the lower the value of g c . Lower 
temperatures favour the insulating phase. 




(12) 



i 4 G 3 n - 2r ] t 2 G 2 n + {r, 2 + P - g 2 )G n - [ V + g{\ - 2w )] = 



(13) 
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The gap opening at g c cannot be directly linked with a characteristic change of the 
ordering of the adiabatic phonons, at least not in the regime of strong coupling and/or 
high temperatures. The anharmonicity plays a crucial role here because it can open a 
gap of width 2g already in the absence of tunneling. This can be compared with the 
situation of harmonic phonons: in equation (0) we replace the distribution Pq by the 
harmonic distribution in (J2J) and study the strong-coupling regime. The treatment of 
the determinant of equation (J5J) in strong coupling does not depend on the type of the 
random variable x r . Thus we obtain the distribution for x r as given in equation (JHJ), 
except for an extra factor P/ l (x r ). In contrast to the distribution of the Ising spins, 
this distribution has only disordered or homogeneously ordered phases for x r but no 
staggered order, implying the absence of a gap. It remains an open question whether or 
not the gap opening and the order of the adiabatic anharmonic phonons are related in 
general. This can be compared with the metal-insulator transition of the Hubbard model 
inside a paramagnetic phase [TUl [TT] . It might be that the gap opening in both cases 
does not require any additional change of order in the models. A common feature of 
these models is the appearence of two eigenstates which are separated by the interaction 
energy if electronic tunneling is absent. This is fundamental also for the gap formation in 
the system with electronic tunneling. Such a separation of electronic state does not exist 
for adiabatic harmonic phonons. A remaining puzzle is that the strong-coupling result in 
equation (jHJ is an antiferromagnetic Ising model that provides a metallic homogeneous 
or a gapped alternating phase, whereas the DMFT result indicates only a gap in the 
disordered phase. 

Experimental relevance of our results can be linked to the complex structures 
observed in manganites and related materials jTHJ ^] . It is believed that the interesting 
properties of these materials (like the colossal magnetoresistance, or CMR) are caused 
by a frozen mixture of ferromagnetic metallic and antiferromagnetic (or charge-ordered) 
insulating clusters. The typical size of these clusters is of the order of ten to hundred 
nanometers. This has been understood as electronic phase separation due to the 
proximity to a first order phase transition. Although our model contains much less 
degrees of freedom than a realistic description of manganites, the phenomenon of self- 
organized phase separation due to frozen metallic (i.e. homogeneously ordered) and 
insulating (i.e. staggered) clusters and a first order phase transition between the 
homogeneously ordered and the staggered phase is also essential in our theory in terms 
of the correlated binary alloy. 

7. Conclusions 

Anharmonic phonons were studied in the adiabatic limit of the Holstein model. We 
applied two different approaches to study the effect of the electron-phonon coupling, one 
is based on a correlated binary alloy (represented by Ising spins), and the other one on 
the dynamical mean-field theory. These approaches cover different regions of the phase 
diagram. In terms of the phonons the correlated binary alloy reveals a homogeneous, an 
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homogeneous phonon phase 



-0.75 




disordered phonon phase 



0.4 0.6 

g 



Figure 1. Mean-field T = 0-phase diagram. Thermal fluctuations are likely to create 
a disordered phase for small values of g and any value of /i. 



alternating (or charge- density wave), and a disordered phase. In terms of the fermions 
there is a metallic state in the homogeneous phonon phase, and insulating states in the 
alternating phonon phase (with a gap) and in the disordered phonon phase (without 
gap). The dynamical mean-field theory, on the other hand, indicates a gapped insulating 
state inside the disordered phase. 
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Appendix A. Effective Ising Models 



The determinant in equation (JHJ) 



P = — det 

Zj 



1 + c 



0(v+t-AgS) 



(A.1) 



is expanded in powers of g (weak-coupling regime) or in powers of the tunneling rate t 
(strong-coupling regime). 
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Figure 2. Density of states for varying coupling stcngths (g = 0.3, 0.4, 0.5, 0.6 and 
1.0) for T = 1, /i = and i = 0.5. Increasing electron-phonon-coupling favors the gap 
formation. 
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0.55 



Figure 3. Chemical potential versus the critical value of the electron-phonon coupling 
g c plotted for t = 0.5 and serveral temperatures T = 0.05, 0.2 and 1.0 (from left to 
right). The lines arc guides to the eyes. 
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In the weak-coupling regime we can expand the argument of the exponential 
function in powers of g. This gives in leading order 



P « exp {-Pg Tr [G e Pifl+t) S] } /N w = exp { - Pg^ [GqS^ 



sMn, 



with Go = (1 + e /3 ^ + ^) 1 and the normalization N w . [Gqq 
Fourier transformation (t — > e) as 



can be evaluated by 



1 + e 



-/3( M +e(fc)) 



-l d d k 
(2tt)° 



In the limit /3 — > oo the integrand becomes the Heavyside step function 0(/x + e(k)). 

The expansion of the expression in equation (jA.l|) in powers of t can be applied in 
the strong-coupling regime: 



exp |Tr [in (l + e 



} « det (G7 1 ) exp 



Tr(GiD) — - Ti{G\DG\D) 



with G^ 1 = 1 + e^- 95) and D = ^ {+9S) - e 
A lenghty but straightforward calculation gives with A r = fi — gS r the relation 



Since there are only values S r = —1,1, the term 



A r -A r , (l + ^)(l + e^A-0 
can also be expressed as a quadratic form with respect to the Ising spins: 

E(S r , S r r) = E + Ei[S r + S r ') + E 2 S r S r ' 

with 



E = - [£7(1, 1) + E{-\, -1) + 2£7(1, -1)] 



(1 _|_ e /3^-/3g)2 _|_ e P/J,+Pgy 



+ 







% 1 + e 2 ^ + 2e^ cosli^g) 



£ 2 = - [E{1, 1) + £(-1, -1) - 2E(l, -1)] 



ft' 



J3n+I3g 



(1 + e^-^9) 2 (1 + e^fa)' 



ft 



e w sinh(/3g) 



and 

Si = l) - £(-1, -l)] = t rr d r , r ^ 

4 o 



Agl + e 2 ^ + 2e^ cosh(pg) j ' 



(1 _|_ e /3M-/3ff)2 (1 _|_ e /3^+/3s)i 
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